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Abstract. This paper has introduced a novel approach to simulate the thermal and 
moisture diffusivity in a semiconducting solid sphere, based on two-temperature theory 
of thermoelasticity. By incorporating two-temperature theory, it accounts for the 
interplay between the temperature and stress, providing a more accurate representation 
of the system’s behavior. The research focuses on the behavior of the sphere when it is 
subjected to a laser pulse that induces varying heat flux on its boundary surface. By 
utilizing the Laplace transform technique, the mathematical model is solved in the 
transformed domain to meet the intended objective. The mathematical model is 
numerically inverted to obtain a comprehensive understanding of the physical 
parameters in the physical domain. A graphic representation of various parameters 
under the effect of moisture diffusivity at two different temperatures are generated using 
the MATLAB software. The model offers a comprehensive approach to accurately 
represent thermal and moisture diffusivities of the solid sphere, enabling a deeper 
comprehension of the various phenomena observed in a wide range of the 
semiconductor devices. By utilizing this model, the engineers and designers can 
enhance their ability to effectively design these devices, ensuring optimal performance 
and functionality. 
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Nomenclature 
bij Kronecker delta, Bij Thermal elastic coupling tensor, 
Cij Strain tensors(mm~!), Ui Components of displacement(m), 
Dm Moisture diffusivity coefficient, aij Two temperature parameter, 
Dg Carrier diffusion coefficients, C€kk Cubical dilatation, 
Ki; Coefficient of Thermal conductivity, p Medium density (Kgm~?), 
€ijk Permutation symbol, DPF, Dr Cross-coupled diffusivities 
A,u ~Lame’s elastic constants, y Conductive temperature, 
Dy Temperature diffusivity K Coupling parameter for thermal activation, 
K ij Materialistic constant, t Time, 
m Moisture concentration Km Moisture diffusivity 
N Carrier density, dy, Coefficient of electronic deformation, 
70 Thermal relaxation parameter, No Carrier concentration atequilibrium position 
On Electronic deformation coefficient, Dg Carrier diffusion coefficient 
at Linear thermal expansion coefficient, Cij Stress tensors(Nm~?), 
Sy surface recombination velocity Cr Specific heat at constant strain, 
c Heat capacity T Photo-generated carrier lifetime, 
mo ‘Two temperatures, To Reference temperature s.t. |T/To| << 1 
7 Thermodynamic temperature, Eg Energy gap of the semiconductor parameter. 


1. Introduction 


The hygro-thermoelastic behavior of a solid semiconductor sphere with two 
temperature is a complex phenomenon that involves the interaction of moisture, heat, 
and mechanical stress within the material. This type of analysis is crucial in 
understanding the performance and reliability of semiconductor devices, which are 
used in a wide range of technological applications, including electronics, solar cells, 
and optoelectronics. The behavior of such semiconductors under different 
environmental conditions, such as varying temperatures and humidity levels, can have 
a significant impact on their functionality and durability. Understanding the hygro- 
thermoelastic behavior of solid semiconductor spheres is crucial for the design and 
optimization of semiconductor devices. By analyzing the effects of temperature, 
humidity, and mechanical stress on the material, we can develop more reliable and 
efficient devices that can withstand different operating conditions. Furthermore, this 
knowledge can help in the development of new materials with improved hygro- 
thermoelastic properties, leading to advancements in various technological fields. 

The study of hygro-thermoelastic behavior in solid spheres under the influence of 
moisture and thermal diffusivity is an important area of research in materials science 
and engineering. Such studies aim to understand the behavior and response of materials 
to coupled effects of temperature, moisture, and mechanical loading. The hygro- 
thermoelastic behavior of a solid sphere can be influenced by several factors, including 
the material properties, moisture diffusion, thermal diffusivity, and boundary 
conditions. 

Moisture diffusion is often modeled using Fick’s law of diffusion, which describes how 
moisture moves within the material due to differences in moisture concentration. The 
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interaction between moisture diffusion and mechanical loading can result in stresses 
and strains, leading to changes in the mechanical response of the sphere. Thermal 
diffusivity describes how heat propagates within a material. In the context of hygro- 
thermoelastic behavior, thermal diffusivity plays a crucial role in determining how 
heat is transferred and distributed within the solid sphere under the influence of 
moisture. It influences the rate at which temperature changes occur and interacts with 
mechanical loading and moisture diffusion. The material properties of the solid 
sphere, such as elastic modulus, Poisson’s ratio, and coefficient of thermal expansion, 
play a significant role in its hygro-thermoelastic behavior. These properties can vary 
with moisture content, which affects the mechanical and thermal response of the 
sphere. Moreover, the boundary conditions applied to the solid sphere can influence 
its response to hygro-thermoelastic effects. For instance, temperature and moisture 
gradients at the boundaries can induce stresses and strains within the sphere, leading 
to deformation or failure. 

Understanding the hygro-thermoelastic behavior of solid sphere under the influence 
of moisture and thermal diffusivity is important for various applications. It can help 
engineers and researchers design and optimize materials and structures that are 
subjected to environmental conditions involving moisture and temperature changes, 
such as composite materials, and electronic devices. 

Researchers in the field of nanotechnology have shown a significant interest in 
semiconductor nanostructures. To fully comprehend — semiconducting 
micro/nanodevices, it is essential to study thermoelasticity. This is because, 
according to thermoelasticity, semiconducting materials can be classified as elastic 
materials. Several theoretical models have been explored to understand the 
relationship between photothermal equations and thermoelasticity. By delving into 
these models, researchers can gain a better understanding of the behavior of 
semiconducting materials and their potential applications in nanotechnology. The 
classical uncoupled thermoelasticity theory, initially proposed by Duhamel [1], 
presents two notable limitations. Firstly, this theory assumes that the state of elastic 
materials remains unaffected by temperature variations, disregarding any potential 
thermal effects on the material’s behavior. Secondly, due to the parabolic nature of 
the heat equation, this theory predicts that temperature propagates at an infinite speed, 
contradicting experimental evidence and the fundamental principles of physics. 
These limitations highlight the need for further advancements in thermoelasticity 
theories to accurately capture the intricate coupling between temperature and 
elasticity in materials. Coupled thermoelasticity, as proposed by Biot[2], offers a 
viable solution to the problems at hand. This theory establishes a relationship between 
the equations of heat conduction and elasticity equations, thereby providing a 
comprehensive understanding of the interplay between these two phenomena. By 
considering the effects of thermal expansion and contraction on the mechanical 
behavior of materials, coupled thermoelasticity enables us to accurately predict the 
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response of structures to thermal loads. This approach has proven to be particularly 
useful in the design of high-temperature ap plications, where the effects of thermal 
stresses can have a significant impact on the performance and durability of materials. 
The theory in question has a limitation in that it only predicts heat waves that 
propagate at an unlimited speed. 

By utilizing Fourier’s law, as suggested by Szekeres [3]-[5], one can effectively 
express the heat flux in conduction as 


Po Bey (1.1) 
The moisture flux potential can be defined by utilizing Fick’s law as: 
ffr= —Dmm,. (1.2) 


Numerous experiments have consistently demonstrated a strong and noteworthy 
relationship between the movement of moisture and the temperature distribution, and 
conversely, the temperature field influences the transport of moisture. It has been 
observed that the existence of moisture concentration gradients leads to the 
emergence of heat fluxes, which are commonly known as Dufour fluxes. These fluxes 
can be precisely defined as follows [3]-[5]: 


q=(q" +9”) = -Dre, — Dim, (1.3) 


Variations in temperature gradients can potentially influence the moisture flux, also 
known as the Soret flux 


f = (f* + fi") = —Dnm, — D9 j- (1.4) 


The potential for heat and moisture flux can be determined by analyzing the diffusion of 
heat and moisture together as 


q=(q" +9?) = —Dry, — Dimi, (1.5) 
f=(f" +f) =-Dnm, — DP p,. (1.6) 


However, when a temperature gradient is suddenly imposed on a homogeneous and 
isotropic medium, Cattaneo [6] and Vernotte [7, 8] propose a more comprehensive form 
of Fourier law that takes into account a relaxation time. In the study conducted by Lord 
and Shulman [9], a comprehensive theory of thermoelasticity was introduced, 
specifically focusing on an isotropic body with a single relaxation time. This theory 
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encompasses a hyperbolic heat equation, indicating that the propagation of temperature 
occurs at a finite speed. The findings of this research shed light on the intricate dynamics 
of thermoelastic behavior, providing valuable insights into the thermal properties of 
materials. Green and Lindsay [10] presented a comprehensive analysis of 
thermoelasticity, offering a more refined comprehension of this field. Their study focused 
on elucidating the symmetry of the linear heat conduction tensor, providing valuable 
insights into the behavior of heat transfer in elastic materials. Dhaliwal and Sherief [11] 
further expanded on this concept by providing comprehensive thermoelasticity equations 
for an anisotropic medium. These findings have significant implications for the study of 
heat transfer and material properties, and can be applied in a variety of fields such as 
engineering and physics. Green and Naghdi [12]-[14] introduced thermoelastic theories 
that encompass both linear and nonlinear behavior, considering the presence or absence 
of energy dissipation. They further extended the Fourier law to account for these 
variations. 

The study conducted by Fernandez & Quintanilla [15] delved into the linear 
thermoelastic deformations of dielectrics. In the scenario where a SM is exposed to 
an external laser beam, the absorption of optical energy leads to the generation of 
excited free electrons with a semiconductor gap energy E g, resulting in the formation 
of a carrier-free charge density. Consequently, both electronic distortion and elastic 
vibration undergo changes. It is important to note that in this particular situation, the 
presence of thermal-elastic-plasma waves will have a significant influence on the heat 
conductivity equations. In their study, Kaur et al. [16] investigated the impact of 
exponential laser pulse on a semi-conducting solid sphere. They specifically focused 
on the influence of MGTPT and Hall current effect on the sphere. In addition to these 
Kaur & Singh [17, 18], Zhang and Li [19], Lata and Kaur [18], Mohamed et al. [20], 
Kaur et al. [21]-[23], Ebrahimi et al. [24], Lata and Kaur [25], Singh et al. [26] also 
done studies on the Hall current effect and other theories of thermoelasticity. 

In this study, we investigate the hygro-thermoelastic behavior of a solid 
semiconductor sphere using two temperature model. This model takes into account 
the coupling between the electronic and lattice temperatures, as well as the effects of 
moisture diffusion and mechanical deformation. By considering these factors, we can 
accurately predict the changes in the material’s properties, such as thermal expansion, 
electrical conductivity, and mechanical strength. In conclusion, the hygro- 
thermoelastic behavior of solid semiconductor spheres with two temperature is a 
complex phenomenon that requires a comprehensive analysis. By considering the 
coupling between temperature, moisture diffusion, and mechanical deformation, we 
can gain valuable insights into the properties and performance of semiconductor 
materials. This knowledge is essential for the design, optimization, and advancement 
of semiconductor devices in various technological applications. 
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2. Basic ecuations 


The equations governing the interaction of hygro-photo-thermoelasticity [3]-[5] are 
given by: 


1. Constitutive relations 


a = (Auk, x _ eT _ B™m _ inN) bij + pb (ui; + Uj,i) (2.1) 
= (3A + 2u)ar (2.2) 
as (3A + 2am (2.3) 
On = (8A + 2p1)dp (2.4) 
y—-T=aV’y (2.5) 
2. Equation of motion 
0ijj + Fi = pi (2.6) 
3. Plasma diffusion equation 
ON N 
= (Deven - ~) + KT, (2.7) 
Ot \ € 
T ONo 
here K = ——— 2.8 
where ~ Or’ (2.8) 
4. Moisture diffusion equation 
Eg N om BF moDm ei; 
[Pin (mj) | + Diy. ij + a | OE K,, (2.9) 
5. Heat conduction equation 
E,N BL To 
D ; | pm P g — 4] 3 2.10 
| T (P35) 4 + Drm, oY On pC” ( ) 


where, DzpCr = Ki = Kidj , i is not summed. The use of a comma after a subscript 
indicates a partial derivative with respect to space coordinates, while a dot above a 
variable denotes differentiation with respect to the time variable f. 
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3. Mathematical model 


The semiconductor solid sphere being discussed (Figure 1) exhibits characteristics 
of axisymmetric, ensuring that it possesses symmetry around its central axis. Ad- 
ditionally, the sphere demonstrates thermal homogeneity, and a radius of 7 . An 
external laser pulse heating device has been utilized to effectively heat the outer 
surface of the sphere. In order to solve the problem at hand, we must take into 
account the given spherical polar coordinate system (7, 6, yg) with the origin at the 
center of the sphere and the z-axis as the axis of symmetry. The initial temperatures 
of the sphere, represented by 7o , have been kept constant and uniform. 


bedded ill te 


Figure 1: Structure of problem 


It is worth noting that, as per the regularity criterion, all fields within a medium are 
assumed to have finite characteristics. This is an important consideration when 
examining the functions in question, as their dependence on both radial distance r 
and time t must be taken into account. This is particularly relevant in the perspective 
of a 1-D problem where symmetry plays a significant role. Therefore, it is crucial 
to carefully consider these factors when analyzing the problem at hand. The 
formulation of the expressions demonstrating the relationships between 
displacement-strain and components of displacement can be achieved as a result of 
the given conditions are: 


U = (Up, Ug, Us) = (u,0,0)(r, t), (3.1) 


Ou u 
——— or’ €00 = Coe = re €ro = Cro = €64 = O. (3.2) 


The equation (3.1) will take the form when utilizing (2.1) as 
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Ore = (A+ 2p) + 2r= ~ (ATT + B™m+ dnN), (3.3) 


Ou 
000 = Cb = \S + 2(A + 1)- _ (°T + B"™m + in) ' (3.4) 


The dilatation term e can be determined by 


—_ 1 A(r*u) (3.5) 


r2. Or 


The equation of motion can be transformed as 
Oru 


Oo . 1 _ 
+ = (20rr — 909 — 949) = Pa: (3.6) 


Or 


Utilizing equations (3.3)-(3.4), in equations (3.6) yields: 


a) 


(A + 2u) 5 


e of O 9 mom ON | (<*) 
~~ Bs (1-eV") p-B"2-- bn =plaz). = B.7) 


The Laplacian operator V7 in the spherical coordinate system, is given by 


ao 20 #120 6 
2). - _ a2 
. = Or? * rOr r20r (: 9) ; (3.8) 


2 0 
Pre-operating both sides of eq (3.7) by € - x) yields 


Or 
i o... gr 2 v2 . Bmyy2 c V2N _ Oe 
(A + 2) Ve — 8° V* (1—aV") y-| m— dn =|(52)- (3.9) 
Also, equations (3.7), (3.9), and (3.10) can be written as, 
ON N 
ON (Ds (v2N) -7) +k(1—aV?) g, (3.10) 
t if 
- ‘ . ‘ E,N . Om - Oe 
Km [DmV?m] + KmnDLV'*¢ ic = = [Kn + B moDn | ; (3.11) 


_ E,N 0 Tp Oe 
KV?9 + pCeD7'V?m + - = [pees (1 av?) e+ "To. (3.12) 


Semiconductor Thermoelastic Solid Sphere under Moisture and Thermal Diffusivity 35 


1 @ 
Pre-operating both sides of eq (3.9) by € F =) we obtain 
2 Ty2 2 my2 cw Pe 
(A + 2u)V*%e — 8° V* (1—aV*) p— B"™ Vm — OnVEN = ap): (3.13) 


The equations mentioned above can be converted into dimensionless form using 
the following subsequent dimensionless quantities. 


l a 
(r’,u’) =von(r,u), (T’,N’, i,,9,m') =—> (87 T, OnN, 015, Bm) , 
pv , 
’ a pCr ‘ in 
(1, 7',t’) = vgn(70,7,t), a = ae n=) pue=rA+2p, y= X42p° (3.14) 


After applying the equation (3.14) to equations (3.10)-(3.13) and removing the 
primes, the resulting equations are obtained as 


Ore 
V7e — V? (1-aV?) p— V?m — (V?N) = (5) ; (3.15) 
ae . 
a (51(V°N) — d2N) + 63 (1-aV ) 9, (3.16) 
’ pen . ON . Om. Oe 
[Vm] + 6504p + 62% = bu mm by | , (3.17) 
2 . 2 ‘ eS 9 . Oe 
V°e + b12V°m + O5N = 0135) (1-aV*) p+ 6a, ; (3.18) 
where 
: ..t «eR = BT E, _ BT? T  . gmpr 
01 = Den, 02 = ’ 03 = T? 05 = z ; 06 = Sk 8 — = ue ; 
T B OnCRr(A + 2p)nT pCRr(A + 2m) vonDm 
E,p™ : 1 .  (6™)?m+0 . pfDm ., pt 
= fone a 010 = ‘ Ou = 5 = ‘ 012 = 53 013 = P 
TV6N7OnKmDm nDm prvponkm vpn? (A + 2p) 


After applying the equation (3.14) to equations (3.3)-(3.4) and removing the 
primes, the resulting equations are obtained as 
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2 Ou 
, Or 


On = +(1-7")e ((1 aV?)p+m+N), 


090 = 84 = fe +(1-7)e-((1-aV?)p+m+N). 


Initial conditions are assumed as 


Ou 
0) =0= —(r,0 
u(r,0) =0 = =*(r,0) 
T(r,0) =0= or We. 0) 
Or 
ON 
N(r,0) =0= Or (r, 0) 
nen=o=— bn) 
Or 


The Laplace transform of a function f w.r.t. time variable t, is defined as 


L(f(t)) = Fls) = | flt)e~tdt. 


) 


Using Laplace transforms from eq. (3.25) to (3.15)-(3.18) yields 


(V? — s?)é-V? (1-aV?) 6- V?m- V?N =0, 
(5, V? — d14)N + 63 (1 -aV?) 6 =0, 
digé + (V? + 629) m + dgV?varphi + dgsN = 0, 
bi5€ + (d17V? + d16) varphi + 65N + 612V?m = 0, 


where 


(3.19) 


(3.20) 


(3.21) 


(3.22) 


(3.23) 


(3.24) 


(3.25) 


(3.26) 


(3.27) 


(3.28) 


(3.29) 


614 = 62 +8, O15 = —dg8, 516 = —6138, 617 = (1+ 46138), b1g = —d118, 529 = —d108. 


Utilizing equation (3.25) to equations (3.19)-(3.20), yields 
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or = Poe + (1-8 ((1—a?V?)6+m+N), (3.30) 
Tm = P= + (1-7?)@-(1-@V*)p+m+N). (3.31) 


When Eqs. (3.26) to (3.29) are decoupled, we get 


(V8 + BV® + CV4 + DV? + E) (é,¢,m, N), (3.32) 


where 


A = —616g512 + 61617 + 61518512 — a61415, 
B = (a63(65 — 859512 — 618612 + 615) + 61458012 — 614517 + 68612514 — 615514 
— 51(—616 — 617629 + (—5g512 + 517)8” — 615 + 618612 + 0515520 — 518617 + 51568))/A, 
C = (—63(65 — 859512 — 518512 + 515) + a53(62955 — 8755 + 512598° — 595158 + 51855 — 515520) 
+ 514(518612 — 616 — 617520 — (51258 — 517)8” — 515 + 2615629 — 618517 + 51558) 
— 51(—616620 — s?(—616 — 517620) — 615620))/A 
D = (—ad362955s8” — 53(52055 — 8755 + 512698° — 596158 + 51855 + 515520) 
+ 614(—616529 — s?(—616 — 617520) — 615520) — 615166208") /A, 
E = (6145165208" + 53520558) /A. 


Let Ai, i= 1, 2, 3, 4 be the roots of equation (2.32), so we can write 


(V? — Az)(V? — AB)(V? — AZ)(V? — AG)(E, F, m, N) = 0 (3.33) 


where \?, i = 1,2,3,4 are the roots of the characteristic equation of equation (3.32) 


dB + B® + Cr + Di? + E=0. (3.34) 


The general solution of Eq. (2.33) can be written in the form 


4 
(é, DP; m, N) = >. Ge Gi Nis Wi) gilty2 (Air) ’ (3.35) 


i=1 


where /,() specifies the modified Bessel functions of second types of order n. The 
relation for Bessel function is 
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7 29/2 fy p(2)de = 2°? Iy/9(2). (3.36) 
By introducing Eq. (3.35) into Eqs. (3.26)-(3.29), we obtain the following expression 


- (Ai — 8”) [(A? + b20)ds — dod128A?] + A?(55d1gs — 515698) — A? [d13612? — 615 (519A? + d20)| 
; (63(1 - ad?) [55.? + 5609 — bg 28)? + (5,4? = 644)s? [53512 = (7 + 599)(517? + 516) , 


(3.37) 


(A2 — s?) [5553 (1 — ad?) — (517d? + 16) (1A? — 514)] — dis (1 — A2) A? [= (610? — 6x4) + (3.38) 


a (63(1 - a?) [65A? + 65629 — 695128A?] + (dj ? = 614) 8? (685124 - (A? + 529) (517A? + 61 
ne (d3(1 — ad—”)) [(A? — 8) (AF + b20 + 5187] (3.39) 
i= 


(53(1 - ad?) [55A? + 65609 — 5g5128,?| + (5143 - 614)8? (ds512\ - (A? + 620) (617A? + 516)) 


The expression for the displacement u in the Laplace transforms domain is as 
follows: 


4 
t= > <-9ilsalXir): (3.40) 


The modified Bessel J, is defined by the following relationships for any positive 


number x. 
f;2.., 
hj = _, sinha, (3.41) 


/ i 3.42 
[372 = = (cosh - a) . ( ) 


By utilizing Eq. (3.41) into Eqs. (3.35), yields 


4 
= [2 ‘ 
(2,9, m, N) = a > 1G me) Sex Sinh(ir). (3.43) 
i=1 . 


The displacement 7 can be expressed in the Laplace transform domain by utilizing 
equation (3.25) in (3.40) as 


= Der Gi ; : sinh(\;r) 
i= V2 7 (cosn(air) — Gur) ) . (3.44) 
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Differentiating equation (3.40) in terms of r yields 


_ 4 
= — "ES gi{l; sinh(A;r) — n; cosh(,;r)}, (3.45) 
ir T 
i=1 
_(2+Ar?) 2 
ae 3 Le Bee S UL 


Thus, using (3.43)-(3.45) in equations (3.30) and (3.31), the expressions for thermal 
stresses are derived as 


2 - ; ~ e 
on = V2 date sinh(A;r) — n1; cosh(A;r)}, (3.46) 


a (3.47) 
500 = \/- i{p; cosh(Ayr) — m; sinh(A;r) }, 
00 Ved {p; cosh(A;r) — m; sinh(A;r) } 
oe (7 1-7? -(G+m) (3.48) 
Pi = 23/2? i 75/22 p1/2y}/? ’ 
— -(G+n 3.49 
l= 7,4 27+) , nl =n. ( ) 
rl/2y1/? 


4. Boundary conditions 


The mechanical boundary condition can be expressed as the exterior surface of the 
sphere being traction-free and constrained by a time-dependent variable heat 


T(R,t) =T,H(t), t > 0, (4.1) 


Ory (R, t) = 0. (4.2) 


In the diffusion phase, carriers have the possibility of reaching the surface of the 
sphere and may recombine with a finite probability. As a result, the boundary 
condition for carrier density is established as 

aN 


De— = 
FOr 


SN, at r=R, (4.3) 


The following moisture diffusion boundary condition is applied at the boundary 
surface when r = ro 
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m=Oatr=ro. (4.4) 


By applying the Laplace transform on (4.1)-(4.4) yields 
To 


(1- aV’) Q(R, s) = = (4.5) 
ary (ts) = 0; (4.6) 
ON 4.7 
E ar = syN(R, s) ( 
m=0 at r=R. (4.8) 
Equations (3.39) and (3.36) are substituted into Eq. (4.5)-(4.8), giving 
Gi To 
V2 xa {( 1— a?) TR sinh(A; R)} = ra (4.9) 
4 
wr { 11; sinh(A;R) — nl; cosh(A;R)} = 0, (4.10) 
i=1 
4 
S— gin { a; cosh(\;R) + 6; sinh(A;R)} = 0, (4.11) 
i=1 
{2d 0h Ni on AiR) = 0, (4.12) 
_ Dpd? —Dr Sy 


» = —— ; 

R1/2 . 2 R3/23/2 AiR 
i 1H ee 2+ d?R? 

ni = hi 7 : = ’ hij = ; 5/2 : 
N20, aoa R3); é 


The values of gi, i = 1, 2, 3, 4 can be obtained by solving equations (4.9)-(4.12) by 
Cramer’s rule 


Aj 
gi(s) = a? (4.13) 


Semiconductor Thermoelastic Solid Sphere under Moisture and Thermal Diffusivity 41 


Gi, Gi2 Giz Gis 
Go, Go2 Go3 Gog 


A= ; 
Gs Geo Gos Gs Gis) 
G4, G42 G43 Gag 
Aj is obtained from A by replacing its i” column by | a 00 0}. (4.15) 


G -/2 Si_(1 — a2) 
i= = VR QA; ; 


Go; = {l1; sinh(A;R) — n1; cosh(A;R)}, 
G3i = ni {a; cosh(A;R) + bj sinh(A;R)}, 


G4i = V7 - sinh(A;f), 


4=1,2 


Using the values of gi(s) from (4.13) in (3.43), (3.44), and (3.46)-(3.47) the 


expression for displacement, carrier density, stresses, and temperature distribution 
can be obtained. 


5. Inversion of the transforms 


To obtain the solution in the physical domain, the transforms in equations (3.49) 
and (4.4) need to be inverted by using 


+ioco _ 
f(z,t)=— - F(a, s)e~*ds. (5.1) 
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To evaluate the integral in equation (5.1), it is recommended to apply Romberg’s 
integration [27] with an adaptive step size. This method involves dividing the 
integration interval into smaller subintervals and approximating the integral using 
a combination of these subintervals. The adaptive step size ensures that the 
approximation is accurate and efficient, as it adjusts the size of the subintervals 
based on the function’s behavior. By using this method, the integral can be 
evaluated with 


6. Numerical results and discutions 


The physical data of silicon (Si) material is utilized in this study to show the 
theoretical findings and visually represent the impact of two temperatures using 
MATLAB software. 


42 Eduard-Marius CRACIUN, Iqbal KAUR, Kulvinder SINGH 


A = 3.64 x 10!°Nm-?, To = 300 K, 

pu = 5.46 x 10!°°Nm-2, Hp = 1Jm='nb-}, 

B=7.04 x x 10® Nm~? deg™!, 7=5x107*s, 

bn = -9 X1072!m-8 , No = 10?°m-3, 

p = 2.33 x 10° Kgm-3, €9 = 8.838 x 10-1!2F' m=}, 

Cre = 695J Kg K-}, Eg = 1.11 e V, 

Kk =150Wm'!K-}, ar =3x 10-°K—}, 

K* = 1.54 x 10?Ws , Sy = 2ms7}, 

Dg = 2.5 x 1073m?2s—! , Ho = 108Col.cm-'s7}, 

un = 4a x 1077Hm"! , on = 9.36 x 10°Col?C—!m-!s—}. 


The data presented in the figures pertains to four distinct cases of the two 
temperatures parameter. These cases are represented by solid black lines (a = 0.02), 
solid red lines (a = 0.04), solid blue lines (a = 0.06), and solid purple lines (a = 
0.08). The analysis of these cases is based on the temperature parameter, which is 
a crucial factor in determining the behavior of the system under consideration. The 
radial displacement u for the hygro- photo-thermoelasticity model with two 
temperatures parameter is illustrated in Figure 2, exhibiting a distinct pattern. 
Initially, the displacement starts from a positive zero value and experiences a rapid 
surge, eventually reaching a peak maximum value near r = 0.7. This behavior can 
be attributed to the combined effects of photo-excitation and moisture diffusivity. 
Notably, as the two temperatures parameter increases, there is a noticeable increase 
in the radial displacement. This observation suggests a correlation between the two 
temperatures parameter and the magnitude of radial displacement, indicating that 
higher values of the two temperatures parameter result in higher value of 
displacement. For the parameter value a = 0.02 of the two temperatures parameter, 
the radial displacement exhibits the lowest deviation. 


Figure 2: The radial displacement variation with two temperatures 
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Figure 3: The variation in temperature distribution with two temperatures 
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Figure 4: The moisture diffusion variation with two temperatures 
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Figure 5: The variation in dilatation distributions with two temperatures 
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Figure 6: The deviation in carrier density with two temperatures 
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Figure 7: The deviation in radial stress with two temperatures 
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Figure 8: The deviation in hoop stress with two temperatures 
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The moisture diffusion behavior of the hygro- photo-thermoelasticity model is 
depicted in Figure 4, considering different values of two temperatures parameter. 
It is noteworthy that a higher value of this parameter leads to increased variability 
in moisture diffusion, while a lower value results in reduced variability. In each of 
the four situations corresponding to the two temperatures parameter, the moisture 
diffusion initiates from zero and attains maximum values near the surface. 
Additionally, it is important to mention that at r=0 and r=R, the moisture diffusion 
complies with the boundary condition specified in equation (2.8) by being zero. 
Figure 5 displays the dilatation distributions for the hygro- photo-thermoelasticity 
model at various values of two temperatures parameter. The results indicate that 
lower values of the two temperatures parameter exhibit greater variability in 
dilatation distributions compared to higher values. Additionally, as the radial 
distance r increases, the variability in dilatation distributions is observed to 
significantly in- crease. These findings suggest that the two temperatures parameter 
plays a crucial role in determining the dilatation distributions of the hygro-photo- 
thermoelasticity model, and that the radial distance r should also be taken into 
account when analyzing these distributions. In each of the four situations 
corresponding to the two temperatures parameter, the dilatation initiates from zero 
and attains maximum values near the surface. The fluctuation in carrier density for 
the hygro-photo-thermoelasticity model at various temperatures is illustrated in 
Figure 6. It is observed that an increase in temperature values leads to an increase 
in carrier density. 

Additionally, it is noted that as the radial distance r expands, there is a significant 
increase in the variance of carrier density. In all four situations corresponding to the 
temperature parameter, the carrier density starts from zero and reaches maximum 
values near the surface. This behavior can be attributed to the combined effects of 
photo-excitation and moisture diffusivity. 

Figures 7 to 8 illustrate the changes in stress components for the hygro-photo- 
thermoelasticity model when considering the two temperatures parameter. It is 
evident that there is an upward trend in the stress components as the values of the 
two temperatures parameter increase. In all four scenarios associated with the two 
temperatures parameter, the o;, starts from zero and reaches its peak values near the 
surface. It is worth mentioning that at r=0 and r=R, the orr adheres to the boundary 
condition specified in equation (2.6) by being zero. Moreover, in each of the four 
situations corresponding to the two temperatures parameter, the 09 begins from 
zero and attains its maximum values near the surface. 


7. Conclusions 


In this study, a unique methodology is employed to explore the influence of 
photothermoelasticity theory on moisture diffusivity in an elastic semiconductor 
sphere. The primary objective of this study is to analyze the effects of applying an 
exponential laser pulse on the surface of an infinite semiconducting solid sphere. 
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The governing equations are derived for hygro- photo-thermoelasticity with two 
temperatures, incorporating both thermal and moisture diffusivity. 

e The study presents graphical illustrations that effectively depict the 
influence of two temperatures parameter on multiple components, including 
displacement, thermal stresses, dilation distributions, moisture diffusivity, 
carrier density, and temperature field. Through these visual representations, 
it becomes evident that the factors associated with the two temperatures 
parameter hold substantial importance in the domains under investigation. 
The graphs provide a comprehensive understanding of the impact of these 
parameters on the various components, thereby contributing to the overall 
comprehension of the subject matter. 

e The thermal effect arises when semiconductor materials are exposed to 
concentrated laser beams or beams of sunlight. These materials find wide- 
ranging applications in the renewable energy sector, with a particular 
emphasis on the solar cell industry, where semiconductor materials play a 
crucial role. 

e The semiconductor sphere with two temperature is an innovative and 
promising technology that has the potential to revolutionize the field of 
thermoelectrics. This concept involves utilizing hyperbolic materials with 
unique energy transfer properties to create a thermal diode effect, enabling 
efficient conversion of heat into electricity. By maintaining a temperature 
gradient across the semiconductor sphere, the device can harvest waste heat 
or solar radiation and convert it into usable electrical power. The two 
temperature approach enhances the efficiency of this process by effectively 
managing the temperature distribution within the material. This technology 
holds great promise for various applications, such as waste heat recovery in 
industrial processes, solar energy harvesting, and even powering electronic 
devices. Through further research and development, the semiconductor 
sphere with two temperature could contribute to a more sustainable and 
clean energy future. 
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